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Abstract 

In this paper we consider a variation of the Merton’s problem with added stochastic volatility 
and finite time horizon. It is known that the corresponding optimal control problem may be 
reduced to a linear parabolic boundary problem under some assumptions on the underlying 
process and the utility function. The resulting parabolic PDE is often quite difficult to solve, 
even when it is linear. The present paper contributes to the pool of explicit solutions for stochastic 
optimal control problems. Our main result is the exact solution for optimal investment in Heston 
model. 


1 Introduction 

The optimal control theory for stochastic processes plays a prominent role in financial mathe¬ 
matics and allows to formulate and solve problems in optimal investment, optimal trading, and 
other areas related to finance. The seminal paper by Robert Merton PQ was generalized in several 
directions. The subsequent research has introduced more realistic and complex asset dynamics 
and more realistic market models [¥]. Similar mathematical problems arise in algorithmic 
trading [5], and market microstructure studies [6]. 

In this paper we consider a variation of the Merton’s problem with stochastic volatility and 
finite time horizon. As shown in [5], the optimal control problem may be reduced to a linear 
parabolic boundary problem under some assumptions on the underlying process and the utility 
function. Generally, due to the lack of solution smoothness for this boundary problem, one needs 
to find a viscosity solution to a non-linear equation. 

The resulting parabolic PDE is often quite difficult to solve, even when it is linear. In several 
special cases the explicit solutions were obtained, see [7] for the model similar to Heston’s, and 
[8] for Cliacko-Viceira model. 

The present paper contributes to the pool of explicit solutions for stochastic optimal control 
problems. Our main result is the exact solution to the optimal control problem within the 
framework of the celebrated Heston model [5]. The method of affine representations, as used 
in 0, cannot be applied here. To obtain the desired explicit solution we employ the theory 
of integral representations and special functions. Quasi-analytical solutions were obtained by 
asymptotic and perturbation methods in ||5] for optimal control of more general processes. 
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The paper is organized as follows: in section [5] we formulate the problem; in section [3] we 
discuss the results of (2| on the representation of Bellman function as the solution of the linear 
parabolic equation, and extend their results to the case of exponential utility; in section [4] we 
present the main results of the paper; in section [5] we analyze the obtained explicit solution; in 
Appendix (section [6]) we provide some technical asymptotics relations. 


2 Formulation of the problem 

Let (X,V) = {X s ,V s ) s >t be a vector stochastic process, given by the triangular system of 
stochastic differential equations 

dX/X = pdt+WdB 1 , X t = x (2.1) 

dV = k(e-V)ds + aWdB 2 , V t = v. 

where B 1 = ( Bj) s > t B 2 = ( B 2 ) s >t are correlated Wiener processes with correlation coefficient 
p, i.e. ( B 1 ,B 2 )t = pt , and p, k, 0, a are some constants. The model for the asset price X 
defined by (12.11) is very well known in financial mathematics, and is called the Heston model [9]. 

In order to solve the system of SDE ’s (EH) one needs to solve the second equation first to 
obtain V, and then solve for X. Process V = (V s ) s >t represents the volatility of asset X. Process 
V is a Feller process, also called a CIR (Cox-Ingersoll-Ross) process in financial mathematics .10) 

. We assume that the Feller condition 2/c0 > <j 2 is satisfied for V. Thus, V is always positive 
given Vo > 0. 

Let us consider a controlled process W = (W s ) s > t given by 

dW = adX , W t = w. (2.2) 

where control a = (a s ) s >t is adapted to filtration T>t = cr{B s — B t , s >t} and J f T a 2 dt < oo. 
Our aim is to find the control on which functional EI7 (Wt) achieves its maximum 

J(w,x,v,t) = supEC/(Wr). (2.3) 

a 

Function J(w,x,v,t) is the Bellman function, which is a solution to HJB ( Hamilton-Jacobi- 
Bellman) equation. Utility function U ( w) can be a power or an exponential utility function. For 
the convenience of notation we will omit the indices. 

w y e~ cw 

Up(w) = —, 7 < 0, Ue(w) = 1-, c > 0. (2.4) 

7 c 

Note, that the logarithmic utility is a special case of the power utility with 7 = 0. This well 
known observation can be checked by taking the limit lim 7 _>.o('a ;7 — l)/7 = log(w) and noting 
that utility functions are defined up to an additive constant. 

The financial interpretation of the processes above is the following: process W represents the 
wealth dynamics of the investor, while control a is the evolution of the investor’s position in the 
asset. The utility function U represents the investor’s preferences. Thus, to invest optimally, 
one should find the optimal investment rule at each time t such that the expected utility of the 
terminal wealth at time T is maximised m • 

Function J(w , x, v, t ) is the solution of Hamilton-Jacobi-Bellman (HJB) equaion. The solution 
of HJB equation in the general case is a non-linear PDE. However, this equation may become 
linear under some choice of utility functions. Zariphopoulou [2] has proposed a method based on 
viscosity solution technique. Under this method the equation (12.31) may be reduced to a linear 
parabolic equation. The method may be applied to systems of stochastic differential equations 
of “triangular type”. Under this restriction the drift and the diffusion of the underlying asset X 
may depend on the volatility process V in arbitrary way. 
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It is much easier to perform qualitative and quantitative analysis of the optimal solution 
while provided with the exact formula. Nevertheless, it is not always possible to find an explicit 
solution to a linear parabolic PDE. Explicit solutions were obtained in some special cases: see [7] 
for the model similar to Heston’s, and [8] for Chacko-Viceira model. Either of the above results 
are based on affine representations of solutions. This technique is not applicable in the model 
we consider (Em¬ 
in the present paper we propose a method based on Laplace transform. This technique allows 
us to obtain an explicit solution in terms of confluent hypergeometric functions. 


3 T. Zariphopoulou’s result and extensions 

A quite general model for optimal investment (12.31) with power utility function Up was considered 
in [2] • The asset price dynamics X dependent on random process V was given by the triangular 
system of stochastic differential equations 


dX/X = fi(V,t)dt + (j(y,t)dB l , 

dV = b(V,t)dt + a(V : t)dB 2 , 


(3.1) 


where B 1 = (Bl) s > t B 2 = (Bf) s > t are correlated Wiener processes with correlation coefficient 
P- 

We assume the necessary restrictions on the coefficients of system ED being satisfied in 
order to guarantee the existence of a strong solution (X, V). 

Theorem 3.1 (Zariphopoulou). For process \3.1\) and power utility function we have the 
following representation of Bellman function i2.3\) 

Jp(w , x , v, t) = —f 1/s (v, t), (3.2) 

7 

where 

6=l+p 2 -^, (3.3) 

1-7 

and function f = (v,t) Is a solution to the following parabolic boundary value problem 


ft + tU 2 (v, t) fw + ( b(v, t)+p 


JP-i'v,t)a(v,t) 
( 1-7 )a(v,t) 


fv + 


7 $ P 2 (v,t) 


f = 0 

1 — 7 2 a 2 (v, t) 

f(v,T) = 1. 


The optimal synthesis is given by 


ctp(w, x, v, t) 


w ( p(v,t) p a(v, t) f v (v,t) \ 
z(l-7) 6 a(v,t) f(v,t) ) 


(3.4) 


(3.5) 


To find the optimal control at time t one should substitute the values for processes W and 
(X, V) into synthesis (13.51) . Note, that representation (13.21) can be explained as follows. Firstly, 
one can separate variables in the Bellman function J(w,x,v,t) = u(w)R(x,v,t). This trick 
is well known for the power utility case. Secondly, from the form of the coefficients and the 
boundary condition, one can see that the solution does not depend on variable x. Thirdly, after 
the above manipulations, the remaining non-linear equation can be linearised by substitution 
R S {v,t) = f(v,t). 

We use similar reasoning as above to extend Zariphopoulou’s formula to the case of exponen¬ 
tial utility. It is our first result in this paper. 
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Proposition 3.1. Under exponential utility function Ue, the Bellman function 12.3 \ ) is repre¬ 
sented by 

e ~ cw 

Je{ x, w,v,t) = 1 - f 1 /s (v, t), (3.6) 

c 

where 

8=1-p\ (3.7) 

and function f = f(v,t ) is a solution to the following parabolic boundary value problem 

Wf = (i8) 

f(v,T) = 1. 


ft + 2 a2 (. v it)f vv + ( b{v,t) - p 


p(v, t)a(v, t) 
a(v,t) 


, 5 p* 
Iv 2 a 2 


The optimal synthesis is given by 


a* E (w,x,v,t) 


J_ / mCM) 

cx \a 2 (v,t) 


P a(v,t) f v (v,t) \ 
8 a (M) f(v,t) ) 


(3.9) 


Note, that constant 5 given by (EZD and equation (13.81) for / can be obtained from power 
utility formula (13.31) and (13.41) by taking limit 7 —> — 00 . 


4 The main result 


Let us return to the original problem. Consider the Heston model given by EH). Using the results 
from section [3l in particular substitutions (13.21) and (13.61) from Theorem 13.II and Proposition 13.II 
respectively, we can reduce the original boundary value problem for Bellman function to the 
following boundary value problem for function / = f(v,t) 

— 77 -fw +(k&- -—-pa - f v - —f + ft = 0. (4.1) 

f(v,T) = 1, 


Equation (14.11) is the same for both cases of power and exponential utility functions, with the 
only difference in the values of constants 6 and C. In the case of power utility function 


6 = 1 + p' 


2 7 


1-7 5 


For the exponential utility case 

8 = 1-p 2 , 

Lemma 4.1. The solution of EZP is given by 


C = -~^0-8. 
1-7 2 


C=^-S. 


/(«, t) = r( r (2 ^+y 2) e ~* M)/2 (*(«. *))■ x (*(«> m, 


where 


T(u,t) 


2 kv 

fj2 — l) 



(1 ~ 8 )p 
pa 



M\,ti{z) is a Whittaker’s function, andF(z) is a Gamma function. 
Moreover, we have 

fv _ (v + A + 1/2) Mi + x,n (^(f) t)) 

f v 


2 C 


(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 
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PROOF. Making substitutions 


- 2fcu , , - 

V = — o-. k{T-t)=T, 

<j z 


f(v,t ) = e~ xf v x e^ 2 h{i 


(4.7) 


we get the following boundary value problem for h 


1 1/4 — rf 


~4 + 


h = —hf . 
v 


h{v, 0) =v- x e-*/ 2 . 

Let G(u;£) be the Laplace transform of function h(v,f) with respect to f 


pOO 

G{v',() = / e < ’ T h(v,T)df. 

Jo 


By denoting \(u) = v 1 A e u / 2 , one can see that 


1 C 1/4 — ?7 2 


4 v 


G = ~x(v)- 


(4.8) 


If made homogeneous, equation (14.81) is a Whittaker equation. Whittaker functions M-^ tV (v) 
and W-£ tV (v) are two linearly independent solutions of a Whittaker equation. One can check by 
substitution that the formula below is the solution of equation (14.81) 


G(v-, 0 = 


r(i /2 + c + ? 7 ) 


L(1 + 2?y) 

To proceed further we use formula 6.669.4 from m 

i 

2 : 


pV poo \ 

M-C,v(v) / x(<P)W- ( + / x(<p)M-c, v (<p)d<p 

JQ J v J 


(4.9) 


T (k+V-v) 


e 2( ai_|_a2 )^ coshx co th 2iy ( -x ) I 2ft (V a i«2sinha;) dx = v ' 2 ' r ' . W v ^(ait)M„^(a 2 t), 

t^a\a 2 i (1 + Zfj,) 


Re ( - + /i — v ) > 0, Re^r > 0, a\> a 2 . 


The above formula from m allows us to rewrite (14.91) as 


G(v\ C) = Vd 


OO poo 


JO JO 
Applying 6.643.2 from [I 


tp 1 / 2 A e 2 V cosh b tanh 2 ^ ( ^ ) I 2v sinh. if} dtpdif. (4.10) 


pOO 

/ x^e~ ax I 2v (2^)dx = 
Jo 


r (n + v + i) , p / « 2 

rV + D r e ° a M --- 


to the internal integral in (14.101) . we obtain 


G(u;C) = e^^^fe^-^tanh- + J j " M A , 

By changing the variable of integration if in (14.121) as 

2 dif 


, Re + ^ >0 

(4.11) 

/ cosh if — 1 _ 
* 

(4.12) 


2 log 


tanh ( $ 


sinh if 


= du, cosh if — 1 = 


1 - e 1 '’ 


2 dif 
sinh if 
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we get 


G (5;C) = e 


-5/2 r(?y — A + 1/2) 

r(i + 27?) 


exp 


2 e 1 ' - 1 


O' 


e" — 1 


Mx, v v 


e v — 1 


ch/. 


Inverting the Laplace transform, and using substitutions (1X71) . we recover the formula for f(v,t) 


/(M) 


i r (77 — A + 1 / 2 ) r N+io ° r° 

2ni r(l + 2rf) J N - ioo J 


v e v 
2 e" - 1 


e CU+U 



\ 


M 


\,T] 


where N is a number such that all residues of the integrand are to the right of it. 
Using the well-known representation of Dirac function 



e" — 1 j 
(4.13) 


1 

2ni 


rN+ioo 

IN —200 


e zC d( = S(z), 


and changing the order of integration in (14.131) . we get 


/(M) = 


r(?7-A + l/2) 
r(l + 2 r?) 


V c 

<5(^ + f)exp<( -- 


2e u -1 


e v — 1 


M\ >Ti 


e v — 1 


dv. 


(4.14) 

Note, that f > 0. Thus, we can complement the range of integration in (14.141) to the whole 
line, and, using the definition of Dirac’s function, namely j(( — z)g{Cf)d(, = g(z) for any 
continuous g, we get 


/(M) 


r(j 7 -A + l/2) f v e f T 
T (1 + 2rj) 6XP { 2 e f - 1 J 




(4.15) 


Finally, by reversing the change of variables (1X71) . we receive the main formula (TOD . The 
expression (14.61) for /,,// is obtained by using differential rules for Whittaker functions (see XT]). 
A 


The following theorem is the analogue of the results from section [3] but for Heston’s model. 
It gives an exact solution for the optimal control and Bellman function. 

Theorem 4.1. For power utility j2.4l) and process \2.1\) the Bellman function 12.31) is given by 


J P (w,x,v,t) = —f 1/5 {v,t), 
7 

and the optimal control is 


(4.16) 


where 

f(v,t) 

6 

^(M) 


a* P {w , x, v, t ) 


w fp pa {r] + X + 1 / 2 ) Mi+\, v ('T'(M)A 
x(l- 7 )Vu S v My, ($(«,t)) / 


(4.17) 


S | A ± iA e -»<«.«)/2 (*(„,()/ Mx,, (*(«,*))■ 


i+^_ 2 _, A = -l?+<iAk, 

1 — 7 < 7 Z pa 

2 kv 


c = -r L - y 6 ' 

1 — 7 2 


a 2 (e k ( T *) — l) 



is a Whittaker’s function, andT(z) is a gamma function. 


dfdv, 
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Theorem 4.2. For exponential utility 1 2.4\) and process \2.1\) The Bellman function is given by 


J E {x , w,v,t) = 1-— f 1/S (v , t .), 


and the optimal control is 


a* E {w , x , v , t) 


J_ //f + pg (?? + A + 1/2) (^(u.^U 

cx \v S v (\H(v,i)) ) 


(4.18) 


(4.19) 


where 


f(v,t ) 
<5 

’I'(M) 


r(p-A + i/2) 
r(2p +1) 


e -« , ( U,t)/2 ( ^f( t ; ) ^ )) A M\ >v (ty(v,t)) , 


1-P 2 



2kv 

(j2 (e fe ( T_t ) — l) 


(1 ~ 
pa 



1 = 



2(7 


M\ tV (z) is a Whittaker’s function, and r( 2 ) is a gamma function. 

Comment. The model (13.11) with drift p(X t ,t) = V t is studied in [7]. For that particular model 
the coefficient for / in (13.41) is proportional to v. This is the reason why the affine representation 
is possible in that specific case. Moreover, if the boundary condition is not constant, the affine 
representation is not possible again. Our method, on the contrary, allows us to obtain the 
required solution for any boundary conditions and for both models, either Heston’s or the model 
considered in [7]. 


5 Analysis 

5.1 Analysis of the optimal control 

Optimal control in (14.171) and (14.191) is proportional to the sum of two terms: the first term p/v 
corresponds to the static portfolio optimization problem. It it simply the ratio of the instanta¬ 
neous drift to the instantaneous variance and does not depend on time or on the parameters of 
the volatility process. The second term 

p a (rj + A + 1/2) M 1+X , v ('Sjv, t)) 

5 v M AlI7 (ttM)) 

represents the hedging of the future opportunity set ( see Merton [14]). This term vanishes in a 
number of important special cases: 

p = 0: no correlation between the two driving Brownian motions, 

7 = 0: when the investor has log-utility. 

The case of zero volatility of volatility limit of Heston model (u=0 in (12.11) 1 turns the asset 
process into a process similar to a geometric Brownian motion but with a deterministically time 
varying volatility. If in addition we set Vo = 0 in (12.11) . then we get a geometric Brownian 
motion. 

Using the asymptotic results from Appendix [S] in case er ~ 0, we get 


f(v,t) ~ 1, 


/»('M) ^ C k(T-t) 

f(v,t) kv 2 


(5.2) 
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Therefore, for a ~ 0 the Bellman function is approximately equal to the utility function, namely 
Jp ~ Up, Je ~ Ue\ and the optimal controls are 


a* P {w, x, v, t) 


a* E (w, x, v, t) 


w f n pa Cp 
x(l — 7 ) \ v Sp kv 2 6 

1 f P | P<r Ce c k(T-t) 
cx \v 5 e kv 2 



where constants Cr_] and <5r _] depend on the choice of utility. 


(5.3) 

(5.4) 


5.2 Hedging interpretation via a bond 

Suppose C is given by (S3 or S3 - With substitution C/v = r , f(v,t) = g{r,t) we can rewrite 

S3 as 


-g rr + hr (to — r) g r -rg + g t = 0, 


9 (r,T) = 1, 


where 


b = 


s/C' h C ' (1 + A) ’ "■ <r 2 (l + A)" 

Using the Feynman-Kac formula we get the following expression for g = g(r, t) 


kC 


m = — - 


g(r,t) = Ee-/"^ s , 


(5.5) 


where process r = ( r s ) s >t is a solution to the following stochastic differential equation 

dr = hr(m — r)dt + br 3 ^ 2 dB 2 , rt = C/vt■ (5-6) 


Equation (15.61) describes the well-known 3/2-model of stochastic interest rates (see [T3]-) In this 
model, if r is interpreted as the short rate, value of a zero coupon bond is given by m , see 
m■ Hence, the ratio f v /f can be expressed as 


fy{v,t ) 



e -/ t T r.d«| n =C/ V 


(5.7) 


Substituting m into formulae for optimal control (14.171) and (14.191) one can get an interpretation 
of the hedging of the future opportunity set via a bond. 


6 Appendix. Technical asymptotics relations 

It follows from (14.41) . that function / = f(v,t) can be written as a function of only 4', namely 
/ = /(HI), where T = is given by (14.51) . This representation allows us to study the 

asymptotics of Bellman function and the optimal control by using the asymptotics of Whittaker 
functions (see, for example, ED- 

Indeed, by using the well known formulae below 

M\, v (z) ~ 2 T)+1/2 (1 + O(z)), z-> 0, M\ r] (z) ~ + ^ e z/2 z~ x , z -> 00, 

T( 1/2-A+ 17) 

we get 

~ r( r ( i+ 2^ ) ^ +1/a(w,t) ^ X{v,t) + 0 (*(^))], *(«>*)-^°> ( 6 - 1 ) 
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and 


f{v,t)~ 1, '&(v,t)->oo. 


( 6 . 2 ) 


For f v {v,t)/f we have 


and 


fy{v,t) 


(77 + A + 1/2) 

V 


^(v,t) -> 0, 


fy(v,t) 

f(v,t) 


2 C 1 

cr 2 c>F(i>,f) ’ 


^(u, t) —> 00 . 


(6.3) 


(6.4) 
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The paper is organized as follows: in section [5] we formulate the problem; in section [3] we 
discuss the results of m on the representation of Bellman function as the solution of the linear 
parabolic equation, and extend their results to the case of exponential utility; in section [4] we 
present the main results of the paper; in section [5] we analyze the obtained explicit solution; in 
Appendix (section O we provide some technical asymptotics relations. 


2 Formulation of the problem 

Let (X,V) = (X s ,V s ) s >t be a vector stochastic process, given by the triangular system of 
stochastic differential equations 

dX/X = pdt+WdB 1 , X t = x (2.1) 

dV = k(e-V)ds + aWdB 2 , V t = v. 

where B 1 = (Bj) s > t B 2 = ( B 2 ) s >t are correlated Wiener processes with correlation coefficient 
p, i.e. ( B 1 ,B 2 )t = pt , and p, k, 0, a are some constants. The model for the asset price X 
defined by (12.11) is very well known in financial mathematics, and is called the Heston model 7]. 

In order to solve the system of SDE ’s (EH) one needs to solve the second equation first to 
obtain V, and then solve for X. Process V = (U s )s>t represents the volatility of asset X. Process 
V is a Feller process, also called a CIR (Cox-Ingersoll-Ross) process in financial mathematics [5] 

. We assume that the Feller condition 2 k& > cr 2 is satisfied for V. Thus, V is always positive 
given V 0 > 0. 

Let us consider a controlled process W = (W a ) s > t given by 

dW = adX , W t = w. (2.2) 

where control a = (a s ) s >t is adapted to filtration T>t = cr{B s — B t , s >t} and J f T a 2 dt < oo. 
Our aim is to find the control on which functional EI7 (Wt) achieves its maximum 

J(w,x,v,t) = supE[/(WV). (2.3) 

OL 

Function J(w,x,v,t) is the Bellman function, which is a solution to HJB ( Hamilton-Jacobi- 
Bellman) equation. Utility function U(w) can be a power or an exponential utility function. For 
the convenience of notation we will omit the indices. 

e -cw 

Up(w) = —, 7 < 0, Ue(w) = 1-, c > 0. (2.4) 

7 c 

Note, that the logarithmic utility is a special case of the power utility with 7 = 0. This well 
known observation can be checked by taking the limit lim 7 _>.o('u ;7 — l)/7 = log(w) and noting 
that utility functions are defined up to an additive constant. 

The financial interpretation of the processes above is the following: process W represents the 
wealth dynamics of the investor, while control a is the evolution of the investor’s position in the 
asset. The utility function U represents the investor’s preferences. Thus, to invest optimally, 
one should find the optimal investment rule at each time t such that the expected utility of the 
terminal wealth at time T is maximised m • 

Function J(w , x, v. t ) is the solution of Hamilton-Jacobi-Bellman (HJB) equaion. The solution 
of HJB equation in the general case is a non-linear PDE. However, this equation may become 
linear under some choice of utility functions. Zariphopoulou [T5] has proposed a method based 
on viscosity solution technique. Under this method the equation (12.31) may be reduced to a linear 
parabolic equation. The method may be applied to systems of stochastic differential equations 
of “triangular type”. Under this restriction the drift and the diffusion of the underlying asset X 
may depend on the volatility process V in arbitrary way. 
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It is much easier to perform qualitative and quantitative analysis of the optimal solution 
while provided with the exact formula. Nevertheless, it is not always possible to find an explicit 
solution to a linear parabolic PDE. Explicit solutions were obtained in some special cases: see [7] 
for the model similar to Heston’s, and [8] for Chacko-Viceira model. Either of the above results 
are based on affine representations of solutions. This technique is not applicable in the model 
we consider (EUl- 

In the present paper we propose a method based on Laplace transform. This technique allows 
us to obtain an explicit solution in terms of confluent hypergeometric functions. 


3 T. Zariphopoulou’s result and extensions 

A quite general model for optimal investment (|2.3[) with power utility function Up was considered 
in m The asset price dynamics X dependent on random process V was given by the triangular 
system of stochastic differential equations 


dX/X = fi(V,t)dt + (j(y,t)dB l , 
dV = b(V,t)dt + a(V,t)dB 2 , 


(3.1) 


where B 1 = (Bl) s > t B 2 = (Bf) s > t are correlated Wiener processes with correlation coefficient 
P- 

We assume the necessary restrictions on the coefficients of system (ED being satisfied in 
order to guarantee the existence of a strong solution (X, V). 

Theorem 3.1 (Zariphopoulou). For process \3.1\) and power utility function we have the 
following representation of Bellman function i2.3\) 

Jp(w , x , v, t) = —f 1/s (v, t), (3.2) 

7 

where 

6=l+p 2 -^, (3.3) 

1-7 

and function f = (v,t) is a solution to the following parabolic boundary value problem 


ft + tU 2 (v, t) fw + ( b(v, t)+p 


,7 p{v,t)a(v,t) 
(1 -j)a(v,t) 


fv + 


7 $ P 2 (v,t) 


f = 0 

1 — 7 2 a 2 (v, t) 

f(v,T) = 1. 


The optimal synthesis is given by 


ctp(w, x, v, t) 


w ( p(v,t) p a(v, t) f v (v,t) \ 
z(l-7) 6 a(v,t) f(v,t) ) 


(3.4) 


(3.5) 


To find the optimal control at time t one should substitute the values for processes W and 
(X, V) into synthesis (13.51) . Note, that representation (13.21) can be explained as follows. Firstly, 
one can separate variables in the Bellman function J(w,x,v,t) = u(w)R(x,v,t). This trick 
is well known for the power utility case. Secondly, from the form of the coefficients and the 
boundary condition, one can see that the solution does not depend on variable x. Thirdly, after 
the above manipulations, the remaining non-linear equation can be linearised by substitution 
R S {v,t) = f(v,t). 

We use similar reasoning as above to extend Zariphopoulou’s formula to the case of exponen¬ 
tial utility. It is our first result in this paper. 
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Proposition 3.1. Under exponential utility function Ue, the Bellman function 12.3 \ ) is repre¬ 
sented by 

e ~ cw 

Je{ x, w,v,t) = 1 - f 1 /s (v, t), (3.6) 

c 

where 

8=1-p\ (3.7) 

and function f = f(v,t ) is a solution to the following parabolic boundary value problem 

Wf = (i8) 

f(v,T) = 1. 


ft + 2 a2 (. v it)f vv + ( b{v,t) - p 


p(v, t)a(v, t) 
a(v,t) 


, 5 p* 
Iv 2 a 2 


The optimal synthesis is given by 


a* E (w,x,v,t) 


J_ / mCM) 

cx \a 2 (v,t) 


P a(v,t) f v (v,t) \ 
8 a (M) f(v,t) ) 


(3.9) 


Note, that constant 5 given by (EZD and equation (13.81) for / can be obtained from power 
utility formula (13.31) and (13.41) by taking limit 7 —> — 00 . 


4 The main result 


Let us return to the original problem. Consider the Heston model given by EH). Using the results 
from section [3l in particular substitutions (13.21) and (13.61) from Theorem 13.II and Proposition 13.II 
respectively, we can reduce the original boundary value problem for Bellman function to the 
following boundary value problem for function / = f(v,t) 

— 77 -fw +(k&- -—-pa - f v - —f + ft = 0. (4.1) 

f(v,T) = 1, 


Equation (14.11) is the same for both cases of power and exponential utility functions, with the 
only difference in the values of constants 6 and C. In the case of power utility function 


6 = 1 + p' 


2 7 


1-7 5 


For the exponential utility case 

8 = 1-p 2 , 

Lemma 4.1. The solution of EZP is given by 


C = -~^0-8. 
1-7 2 


C=^-S. 


/(«, t) = r( r (2 ^+y 2) e ~* M)/2 (*(«. *))■ x (*(«> m, 


where 


T(u,t) 


2 kv 

fj2 — l) 



(1 ~ 8 )p 
pa 



M\,ti{z) is a Whittaker’s function, andF(z) is a Gamma function. 
Moreover, we have 

fv _ (v + A + 1/2) Mi + x,n (^(f) t)) 

f v 


2 C 


(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 
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PROOF. Making substitutions 


- 2fcu , , - 

V = — o-. k{T-t)=T, 

<j z 


f(v,t ) = e~ xf v x e^ 2 h{i 


(4.7) 


we get the following boundary value problem for h 


1 1/4 — rf 


~4 + 


h = —hf . 
v 


h{v, 0) =v- x e-*/ 2 . 

Let G(u;£) be the Laplace transform of function h(v,f) with respect to f 


pOO 

G{v',() = / e < ’ T h(v,T)df. 

Jo 


By denoting \(u) = v 1 A e u / 2 , one can see that 


1 C 1/4 — ?7 2 


4 v 


G = ~x(v)- 


(4.8) 


If made homogeneous, equation (14.81) is a Whittaker equation. Whittaker functions M-^ tV (v) 
and W-£ tV (v) are two linearly independent solutions of a Whittaker equation. One can check by 
substitution that the formula below is the solution of equation (14.81) 


G(v-, 0 = 


r(i /2 + c + ? 7 ) 


L(1 + 2?y) 

To proceed further we use formula 6.669.4 from [Sj 

1 

2 : 


pV poo \ 

M-C,v(v) / x(<P)W- ( , v (<p)d<p + W- Cri (v) / x(<p)M-c, v (<p)d<p 

J 0 J v J 


(4.9) 


T (k+V-v) 


e 2( ai_|_a2 )^ coshx co th 2iy ( -x ) I 2ft (V a i«2sinha;) dx = v ' 2 ' r ' . W v ^(ait)M„^(a 2 t), 

t^a\a 2 i (1 + Zfj,) 


Re ( - + /i — v ) > 0, Re^r > 0, a\> a 2 . 


The above formula from .9] allows us to rewrite m as 

G(u;C) = ^^ J <p- 1/2 - A e-¥- i ^ cosh,/ ’tanh 2C / 2 ^ (y^sinhV>) 

Applying 6.643.2 from [5] 

pOO 

/ x^~^e~ ax I 2 i, (2/3\/x) dx = 

Jo 


dipdif. (4.10) 


r (/x + v + i) . „ / r2 

■W r "“ 


to the internal integral in (14.101) . we obtain 


G(5;C) = e^^^fe^-^tanh- ^osh^ + 1 ^^ 

By changing the variable of integration if in (14.121) as 

2dif 


, Re + ^ >0 

(4.11) 

( cosh if — 1 _ 

\-^— v 

(4.12) 


2 log 


tanh ( t 


sinh if 


= du, cosh if — 1 = 


1 - e v ' 


2 dif 
sinh if 
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we get 


G (5;C) = e 


-5/2 r(?y — A + 1 / 2 ) 

r(i + 27?) 


exp 


2 e 1 ' - 1 


O' 


e" — 1 


Mx, v v 


e v — 1 


ch/. 


Inverting the Laplace transform, and using substitutions (1X71) . we recover the formula for f(v,t) 


/(M) 


i r(77 — A +1/2) r N+io ° r° 
2ni r(l + 2rf) J N - ioo J 


v e v 
2e" - 1 


e CU+U 



\ 


M 


\,T] 


where N is a number such that all residues of the integrand are to the right of it. 
Using the well-known representation of Dirac function 



e" — 1 j 
(4.13) 


1 

2ni 


rN+ioo 

IN —200 


e zC d( = S(z), 


and changing the order of integration in (14.131) . we get 


/(M) = 


r(?7-A + l/2) 
r(l + 2r?) 


V c 

< 5 (^ + f)exp<( -- 


2e u -1 


e v — 1 


M\ >Ti 


e v — 1 


dv. 


(4.14) 

Note, that f > 0. Thus, we can complement the range of integration in (14.141) to the whole 
line, and, using the definition of Dirac’s function, namely j(( — z)g{Cf)d(, = g(z) for any 
continuous g, we get 


/(M) 


r(j 7 -A + l/2) f v e f ) 
T(1 + 2rj) 6XP { 2 e f - 1 J 




(4.15) 


Finally, by reversing the change of variables (ED , we receive the main formula (14.41) . The 
expression (14.61) for f v /f is obtained by using differential rules for Whittaker functions (see fl]). 
A 


The following theorem is the analogue of the results from section [3] but for Heston’s model. 
It gives an exact solution for the optimal control and Bellman function. 

Theorem 4.1. For power utility j2-4l) and process \2.1\) the Bellman function 12.31) is given by 


Jp(w,x,v,t) = —f 1/5 {v,t), 
7 

and the optimal control is 


(4.16) 


where 

/(M) 
S 

^(M) 


a* P {w , x, v, t ) 


w fp pa {r] + X + 1/2) Mi+\, v ('f'(M)A 
x(l- 7 )Vu S v )’ 


(4.17) 


S | A ± iA e -»<«.«)/2 (*(„,()/ Mx,, (*(«,*))■ 


1 — 7 <7 Z pa 

2 kv 


c = -r L - y 6 ' 

1 — 7 2 


a 2 (e k ( T *) — l) 



is a Whittaker’s function, andT(z) is a gamma function. 


dfdv, 
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Theorem 4.2. For exponential utility 1 2.4\) and process \2.1\) The Bellman function is given by 


J E {x , w,v,t) = l -— f 1/S (v , t), 


(4.18) 


and the optimal control is 


where 


a* E {w , x, v, t ) 


J_ /ft + per (r] + \ + 1/2) (#(M)A 

CX \V S v £)) ) 


(4.19) 


f(v,t) 

S 


rfo-A + 1/2) 

r(2?? +1) 


e -«r(»,t)/2(^ (t;jt ))A M ^ 


WM)), 


1 -p. 



2kv 

(j2 (e fe ( T -t) — i) 


(1 ~ <5)ft 
Per 



1 = 



2(7 


is a Whittaker’s function, and r(z) is a gamma function. 

Comment. The model (13.11) with drift p,(X t ,t) = V t is studied in [TT] . For that particular 
model the coefficient for / in (13.41) is proportional to v. This is the reason why the affine repre¬ 
sentation is possible in that specific case. Moreover, if the boundary condition is not constant, 
the affine representation is not possible again. Our method, on the contrary, allows us to obtain 
the required solution for any boundary conditions and for both models, either Heston’s or the 
model considered in HD. 


5 Analysis 

5.1 Analysis of the optimal control 

Optimal control in (14.171) and (14.191) is proportional to the sum of two terms: the first term p/v 
corresponds to the static portfolio optimization problem. It it simply the ratio of the instanta¬ 
neous drift to the instantaneous variance and does not depend on time or on the parameters of 
the volatility process. The second term 

p a (rj + A + 1/2) M 1+X , v (^(v, t)) 

5 v M AlI7 (*M)) 

represents the hedging of the future opportunity set ( see Merton m)- This term vanishes in a 
number of important special cases: 

p = 0: no correlation between the two driving Brownian motions, 

7 = 0: when the investor has log-utility. 

The case of zero volatility of volatility limit of Heston model (u=0 in (12.11) '! turns the asset 
process into a process similar to a geometric Brownian motion but with a deterministically time 
varying volatility. If in addition we set Vo = 0 in (12.11) . then we get a geometric Brownian 
motion. 

Using the asymptotic results from Appendix [S] in case er ~ 0, we get 


f(v,t) ~ 1, 


/»('M) ^ C k(T-t) 

f(v,t) kv 2 


(5.2) 
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Therefore, for a ~ 0 the Bellman function is approximately equal to the utility function, namely 
Jp ~ Up, Je ~ Ue\ and the optimal controls are 


a* P {w, x, v, t) 


a* E (w, x, v, t) 


w f n pa Cp 
x(l — 7) \ v Sp kv 2 6 

1 f P | P<r Ce c k(T-t) 
cx \v 5 e kv 2 



where constants Cr_] and <5r_] depend on the choice of utility. 


(5.3) 

(5.4) 


5.2 Hedging interpretation via a bond 

Suppose C is given by S3 or S3 - With substitution C/v = r , f(v,t) = g{r,t) we can rewrite 

S3 as 


-g rr + hr (to — r) g r -rg + g t = 0, 


9 (r,T) = 1, 


where 


b = 


s/C' h C ' (1 + A) ’ "■ <r 2 (l + A)" 

Using the Feynman-Kac formula we get the following expression for g = g(r, t) 


kC 


m = — - 


g(r,t) = Ee-/"^ s , 


(5.5) 


where process r = ( r s ) s >t is a solution to the following stochastic differential equation 

dr = hr(m — r)dt + br 3 ^ 2 dB 2 , rt = C/vt■ (5-6) 


Equation (15.61) describes the well-known 3/2-model of stochastic interest rates (see 0-) In this 
model, if r is interpreted as the short rate, value of a zero coupon bond is given by S3, see [3j. 
Hence, the ratio /„// can be expressed as 


fv(v,t -) 



e -f t T r s ds \ n=C / v 


(5.7) 


Substituting S3 into formulae for optimal control (14.171) and (14.191) one can get an interpretation 
of the hedging of the future opportunity set via a bond. 


6 Appendix. Technical asymptotics relations 

It follows from (14.41) . that function / = f(v,t) can be written as a function of only 4', namely 
/ = /(iff), where T = is given by (14.51) . This representation allows us to study the 

asymptotics of Bellman function and the optimal control by using the asymptotics of Whittaker 
functions (see, for example, HD- 

Indeed, by using the well known formulae below 

M\, v (z) ~ 2 T)+1/2 (1 + O(z)), z-> 0, M\ r] (z) ~ + ^ e z/2 z~ x , z -> 00, 

T( 1/2 —A+ ?y) 

we get 

~ r( r ( i+ 2 ^ ) ^ +1/a(w,t) [^ A («d) + Q(^(t>,*))], *(«,*)-»• 0, (6.1) 
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and 


f{v,t)~ 1, '&(v,t)->oo. 


( 6 . 2 ) 


For f v {v,t)/f we have 


fy{v,t) 


(77 + A + 1/2) 

V 


^(v,t) -> 0, 


(6.3) 


and 


fy(v,t) 

f(v,t) 


2 C 1 

cr 2 yip(i>,f) ’ 


^(u, t) —> 00. 


(6.4) 


Acknowledgements. The authors are grateful to Yury Kabanov for fruitful discussions and 
useful suggestions. 


References 

[1] M. Abramowitz, I. Stegun (1972) Handbook of Mathematical Functions with Formulas, Graphs 
and Mathematical Tables, Dover, New York. 

[2] M. Boguslavsky, E. Boguslavskaya (2003) Arbitrage under power, RISK, 06, 49-53. 

[3] D. Brigo, F. Mercurio (2006) Interest Rate Models - Theory and Practice: with Smile, Infla¬ 
tion and Credit, Springer Finance. 

[4] G. Cliacko and L. Viceira (2005) Dynamic consumption and portfolio choice with stochastic 
volatility in incomplete markets, Review of Financial Studies, 18, 1369-1402. 

[5] J.C. Cox, J.E. Ingersoll, S.A. Ross (1985) A theory of the term structure of interest rates, 
Econometrica, 53, 2, 385-407. 

[6] J. Fouque, R. Sircar and T. Zariphopoulou (2013) Portfolio optimization and stochastic 
volatility asymptotics, Preprint 

[7] S. Heston (1993) A closed form solution for options with stochastic volatility with applications 
to bond and currency options. Review of Financial Studies, 6, 2, 327-343. 

[8] J. Gatheral, N. Taleb (2006) Volatility Surface: a Practioner Guide, Wiley Finance. 

[9] Gradshteyn I.S., Ryzhik I.M. (1980) Tables of Integrals, Series and Products, Academic Press, 
New York. 

[10] Yu. Kabanov, M. Safarian (2009) Market with Transaction Costs. Mathematical Theory, 
Springer-Verlag. 

[11] H. Kraft (2005) Optimal portfolios and Heston’s stochastic volatility model: an explicit so¬ 
lution for power utility, Quantitative Finance, 5, 303-313. 

[12] R.C. Merton (1971) Optimum consumption and portfolio rules in a continuous-time model, 
Journal of Economic theory, 3, 4, 373-413. 

[13] R.C. Merton (1990) Continuous-Time Finance, Blackwell Publishers. 

[14] A. Obizhaeva, J. Wang (2013) Optimal trading strategy and supply/demand dynamics, Jour¬ 
nal of financial markets, 16, 1-32. 

[15] T. Zariphopoulou (2001) A solution approach to valuation with unhedgeable risks, Finance 
and Stochastics, 5, 61-82. 


9 



